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ABSTRACT 

We discuss effects of fluctuation geometry on primordial nucleosynthesis. For 
the first time we consider condensed cylinder and cylindrical-shell fluctuation 
geometries in addition to condensed spheres and spherical shells. We find that a 
cylindrical shell geometry allows for an appreciably higher baryonic contribution 
to be the closure density (fib^o < 0.2) than that allowed in spherical 
q ■ inhomogeneous or standard homogeneous big bang models. This result, which 

is contrary to some other recent studies, is due to both geometry and recently 
revised estimates of the uncertainties in the observationally inferred primordial 
CLi 1 light-element abundances. We also find that inhomogeneous primordial 

nucleosynthesis in the cylindrical shell geometry can lead to significant Be and 
B production. In particular, a primordial beryllium abundance as high as [Be] 
= 12 + log(Be/H) pa —3 is possible while still satisfying all of the light-element 
abundance constraints. 

Subject headings: cosmology: theory - dark matter - early universe - nuclear 
reactions, nucleosynthesis, abundances 
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1. INTRODUCTION 



The analysis of primordial nucleosynthesis provides valuable limits on cosmological 
and particle physics parameters through a comparison between the predicted and inferred 
primordial abundances of D, 3 He, 4 He, and 7 Li. For standard homogeneous big bang 
nucleosynthesis (HBBN) the predicted primordial abundances of these light-elements are in 
accord with the value inferred from observation provided that baryon-to-photon ratio (= 77) 
is between about 2.5 x 10~ 10 and 6 x 10 -10 . This corresponds to an allowed range for the 



baryon fraction of the universal closure density f2 BBBN ([Walker et al. 19911 ; |Smith, Kawano 



|fc Malaney 1993| ; |(Jopi, Schramm fc Turner 1995) ; |Schramm fc Mathews 1995| ), 



0.04<a HBBr V 



50 ~ 



< 0.08, 



where rj = 6.6 x 10 9 £lb h\ Q . The lower limit on f2ft HBBN arises mainly from the upper 



limit on the deuterium plus 3 He abundance (|Yang et al. 1984 ; |Walkcr et al. 1991 ; Smith 



Kawano, fc Malaney 1993|) , and the upper limit to fib arises from the upper limit on the 4 He 



mass fraction Y p and/or the deuterium abundance D/H > 1.2 x 10 5 ( [Linsky et al. 1993| , 
1995 ). Here, /150 is the Hubble constant in units of 50 km s ~ x Mpc~ x . The fact that this 



range for fib hf is so much greater than the current upper limit to the contribution from 
luminous matter fl b Lum < 0.01 (see however |Jedamzik, Mathews fc Fuller 1995] ) is one of 



the strongest arguments for the existence of baryonic dark matter. 

Over the years HBBN has provided strong support for the standard, hot big bang 
cosmological model as mentioned above. However, as the astronomical data have become 
more precise in recent years, a possible conflict between the predicted abundances of the 
light element isotopes from HBBN and the abundances inferred from observations has been 
suggested ( piive fc Steigman 1995| ; [Steigman 1996a| ; [l'urner et al. 1995] ; [Hata et al. 1996| ; 



see also |Hata et al. 1995| ) . 



There is now a good collection of abundance information on the 4 He mass fraction, Y p , 





Pagel 1993 




Izatov, 



Thuan fc Lipovetsky 1994| ; [Skillman fc Kcnnicutt 19"95| ). In an extensive study based upon 
these observations, the upper limit to r\ from the observed 4 He abundance was found to be 
~ 3.5 x 10~ 10 ( piive fc Steigman 1995| ; |Qlive fc Scully 1996j ) when a systematic error in 
Y p of AY sys = 0.005 is adopted. Recently, it has been recognized that the AY sys may even 



be factor of 2 or 3 larger ( |l'huan,, JNature., Izatov, fchipovctsky 19"9"E| ; |(Jopi, Schramm fc" 
Turner 1995t |Schramm fc Mathews 1995| ; (Sasselov fc Goldwirth 1995|) , making the upper 
limit to t] as large as 7 x 10~ 10 . 



On the other hand, the lower bound to 77 has been derived directly from the upper 
bound to the combined abundances of D and 3 He. This is because it is believed that 
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deuterium is largely converted into 3 He in stars; the lower bound then applies if, as has 
generally been assumed, a significant fraction of 3 He survives stellar processing ( [Walker et 
M- 1991| ). 



However, there is mounting evidence that low mass stars destroy 3 He ( |Wasserburg. 



Boothroyd fc Sackmann al 1995| ; |Charbonnel 19951) , although it is possible that massive 



stars produce 3 He. Therefore, the uncertainties of chemical evolution models render it 
difficult to infer the primordial deuterium and 3 He abundances by using observations of the 
present interstellar medium (ISM) or from the solar meteoritic abundances. Recent data 
and analysis lead to a lower bound of 77 > 3.5 x 10~ 10 on the basis of D and 3 He ( pearborn 



Steigman, fc Tosi 1996] ; |Hata et al. 1996| |Steigman 1996a| ; |Steigman fc Tosi 1995| ), if the 



fraction of 3 He that survives stellar processing in the course of galactic evolution exceeds 
1/4. This poses a potential conflict between the observation (Y p with low AY sys , D) and 
HBBN. 



In this context, possible detections ([Songaila et al. 1994 ; Carswcll, Science., et, al. 1994 



19961 ll'ytler fc Fan 1994j ; |lytler, Fan, fc Burlers 19961 ; Hugers fc Hogan 19964 1996b; 



Wampler et al. 1996| ) of an isotope-shifted Lyman-a absorption line at high redshift [z > 3) 



along the line of sight to quasars are of considerable interest. Quasar absorption systems 
can sample low metallicity gas at early epochs where little destruction of D should have 
occurred. Thus, they should give definitive measurements of the primordial cosmological D 
abundance. A very recent high resolution detection by Rugers & Hogan (1996a) suggests a 
ratio D/H of 

D/H = 1.9 ±0.4 x 10 -4 . (2) 

This result is consistent with the estimates made by Songaila et al. (1994) and Carswell 
et al. (1994), using lower resolution. It is also similar to that found recently in another 
absorption system by Wampler et al. (1996), but it is inconsistent with high resolution 
studies in other systems at high redshift (Tytler, Fan & Buries 1996; Buries & Tytler 1996) 
and with the local observations of D and 3 He in the context of conventional models of 
stellar and Galactic evolution ( Edmunds 1994 ; Gloeckler fc Geiss 1996 ). If the high value of 



D/H is taken to be the primordial abundance, then the consistency between the observation 
and HBBN is recovered and the allowed range of fib inferred from HBBN changes to 
^hbbn h 2^ = g Q24 ± 0.002 QJedamzik, Fuller fc Mathews 19941 ; |Krauss fc Kernan 1994 ; 



Vangioni-Flam fc Casse 1995|) . In this case, particularly if h 50 is greater than ~ 1.5, the big 
bang prediction could be so close to the baryonic density in luminous matter that little or 
no baryonic dark matter is required QPersic fc Salucci 1992| ; |Jedamzik, Mathews fc Fullci 



1995|) . This could be in contradiction with observation, particularly if the recently detected 



microlensing events (Alcock et al. 1993, 1994, 1995abc; [Aubourg et al. 1993| ) are shown to 
be baryonic. This low baryonic density limit would also be contrary to evidence ( |White e~t| 
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al. 1993; White and Fabian 1995 ) that baryons in the form of hot X-ray gas may contribute 



a significant fraction of the closure density 

The observations by Tytler, et al. (1996) and Buries & Tytler (1996) yield a low value 
of D/H. Their average abundance is 

D/H = 2.4 ±0.9 x HT 5 , (3) 

with ±2<t statistical error and ±lcr systematic error. This value is consistent with the 
expectations of local galactic chemical evolution. However, this value would imply an 
HBBN helium abundance of Y p = 0.249 ± 0.003 which is only marginally consistent with 
the observationally inferred Y p even if the high AY sys is adopted. 

With this in mind, it is worthwhile to consider alternative cosmological models. 
One of the most widely investigated possibilities is that of an inhomogeneous density 
distribution at the time of nucleosynthesis. Such studies were initially motivated by 
speculation ([Witten 1984| ; [Applegate fc Hogan 1985|) that a first order quark-hadron phase 
transition (at T ~ 100 MeV) could produce baryon inhomogeneities as baryon number 
was trapped within bubbles of shrinking quark-gluon plasma. In previous calculations 
using the baryon inhomogeneous big bang nucleosynthesis (IBBN) model, it has been 
usually assumed that the geometry of baryon density fluctuations is approximated by 
condensed spheres. Such geometry might be expected to result from a first order QCD 
phase transition in the limit that the surface tension dominated the evolution of shrinking 
bubbles of quark-gluon plasma. However, the surface tension may not be large flKajantie 



Karkkainen fc Rummukainen 1990] , 1991, 1992) during the QCD transition, which could 



lead to a "shell" geometry or the development of dendritic fingers (Freese & Adams 1990). 
Furthermore, such fluctuations might have been produced by a number of other processes 
operating in the early universe (cf. |Malaney fc Mathews 199"3| ), for which other geometries 
may be appropriate, e.g. strings, sheets, etc. Thus, the shapes of any cosmological baryon 
inhomogeneities must be regarded as uncertain. 

The purpose of this paper is, therefore, to explore the sensitivity of the predicted 
elemental abundances in IBBN models to the geometry of the fluctuations. We consider 
here various structures and profiles for the fluctuations in addition to condensed spheres. 
Mathews et al. (1990, 1994, 1996) found that placing the fluctuations in spherical shells 
rather than condensed spheres allowed for lower calculated abundances of 4 He and 7 Li for 
the same fib, and that a condensed spherical geometry is not necessarily the optimum. Here 
we show that a cylindrical geometry also allows for an even higher baryonic contribution to 
the closure density than that allowed by the usually adopted condensed sphere. It appears 
to be a general result that shell geometries allow for a slightly higher baryon density. This 
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we attribute to the fact that, for optimum parameters, shell geometries involve a larger 
surface area to volume ratio and hence more efficient neutron diffusion. 

An important possible consequence of baryon inhomogeneities at the time of 
nucleosynthesis may be the existence of unique nucleosynthetic signatures. Among the 
possible observable signatures of baryon inhomogeneities already pointed out in previous 
works are the high abundances of heavier elements such as beryllium and boron ([Boyd 



fe Kajino 19891 |Kajino fc Boyd 199q ; |Malaney fc Fowler 1989| ; Terasawa & Sato 1990; 



|Kawano et al. 1991| ), intermediate mass elements ([Kajino, Mathews fe Fuller 1990|) , or 



heavy elements ( [Malaney fc Fowler 1988] ; |Applegate, Hogan fc Scherrer 1988| ; [Rauscher 



|et al. 19l?4] ). Such possible signatures are also constrained, however, by the light-element 
abundances. It was found in several previous calculations that the possible abundances of 
synthesized heavier nuclei was quite small (e.g., |Alcock et al. 1990| Terasawa & Sato 1990; 



IRauschcr et al. 1994TJ . We find, however, that substantial production of heavier elements 
may nevertheless be possible in IBBN models with cylindrical geometry. 



BARYON DENSITY INHOMOGENEITIES 



After the initial suggestion (Witten 1985) of QCD motivated baryon inhomogeneities 
it was quickly realized (|Applegate fc Hogan 1985| ; |Applegate, Hogan fc Scherrer 1987[) that 



the abundances of primordial nucleosynthesis could be affected. A number of papers have 
addressed this point ( |Alcock, Fuller, fc Mathews 1987| ; [Applegate, Hogan fc Scherrer 198~7j . 



1988; |Fuller, Mathews, & Alcock 1988; 


Kurki-Suonio et al. 1988 


, 1990; 


Terasawa & Sato 


19894 IN, H, (1990; 


Kurki-Suonio & Matzner 1989 


, 199C; 


Vlathews et al. 1990 


, Mathews, 


Schramm & Meyer 1993 




Mathews, Kajino & Orito 199C 


; Jedamzik, Fuller & Mathews] 



BMt |Jedamzik, Mathews fc Fuller 1995| ; |l'homas et al. 1994| ; |Kauscher et al. 1994| ). Most 



recent studies in which the coupling between the baryon diffusion and nucleosynthesis has 
been properly accounted for (e.g., [Terasawa fc Sato 19894 U @ ; |1990| ; [Kurki-Suonio fc 



Matzner 1989|, |1990 


; Mathews et al. 1990 


, Mathews, Schramm & Meyer 1993; 


Jedamzik, 


Fuller & Mathews 1994 


; Thomas et al. 1994 


) have concluded that the upper limit on Qb h 



is virtually unchanged when compared to the upper limit on h 2 derived from standard 
HBBN. It is also generally believed (e.g. |Vangioni-Flam fc Casse 1995| ) that the same holds 



true if the new high D/H abundance is adopted. 

However, in the previous studies, it was usually assumed that a fluctuation geometry 
of centrally condensed spheres produces the maximal impact on nucleosynthesis. Here we 
emphasize that condensed spheres are not necessarily the optimal nor the most physically 
motivated fluctuation geometry. 
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Several recent lattice QCD calculations ( [Kajantic, Karkkaincn fc llnmmukainen 1990| , 
1991, 1992; [Brower et al. 1992]) indicate that the surface tension of nucleated hadron 
bubbles is relatively low. In this case, after the hadron bubbles have percolated, the 
structure of the regions remaining in the quark phase may not form spherical droplets 
but rather sheets or filaments. We do note that the significant effects on nucleosynthesis 
may require a relatively strong first order phase transition and sufficient surface tension to 
generate an optimum separation distance between baryon fluctuations ([Fuller, Mathews,~S 



Alcock 1988|) . However, even if the surface tension is low, the dynamics of the coalescence 



of hadron droplets may lead to a large separation between regions of shrinking quark-gluon 
plasma. Furthermore, even though lattice QCD has not provided convincing evidence for a 
strongly first order QCD phase transition (e.g., [Fukugita fc Hogan 1991| ), the order of the 
transition must still be considered as uncertain QCottlieb 1991| ; |Petersson 1993| ). It depends 



sensitively upon the number of light quark flavors. The transition is first order for three or 
more light flavors and second order for two. Because the s quark mass is so close to the 
transition temperature, it has been difficult to determine the order of transition. At least 
two recent calculations fllwasaki et al. 1995| ; |Kanaya 1996|) indicate a clear signature of a 



first order transition when realistic u, d, s quark masses are included, but others indicate 
either second order or no phase transition at all. 

In addition to the QCD phase transition, there remain a number of alternative 
mechanisms for generating baryon inhomogeneities prior to the nucleosynthesis epoch 
(cf. [Malaney fc Mathews 1993|) , such as electroweak baryogenesis ( [Fuller et al. 199~4[ ), 
inflation-generated isocurvature fluctuations ( polgov fc Silk 1993| ), and kaon condensation 
( [Nelson 1990|) . Cosmic strings might also induce baryon inhomogeneities through 



electromagnetic ( |Malaney fc Butler 1989| ) or gravitational interactions 



Since the structures, shapes, and origin of any baryon inhomogeneities are uncertain, 
a condensed spherical geometry is not necessarily the most physically motivated choice. 
Indeed, we will show that a condensed spherical geometry is also not necessarily the 
optimum to allow for the highest values for fi& while still satisfying the light-element 
abundance constraints. Here we consider the previously unexplored cylindrical geometry. 
String geometries may naturally result from various baryogenesis scenarios such as 
superconducting axion strings or cosmic strings. Also, the fact that QCD is a string theory 
may predispose QCD-generated fluctuations to string-like geometry ( |Kajino fc Tessie 1993| ; 
lassie fc Kajino 1993| ). Hence, cylindrical fluctuations may be a natural choice. 
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3. OBSERVATIONAL CONSTRAINTS 

We adopt the following constraints on the observed helium mass fraction Y p and 7 Li 
taken from Balbes et al. (1993), Schramm & Mathews (1995), Copi et al. (1995) and Olive 
(1996): 

0.226 <Y P < 0.247, (4) 
0.7 x 10~ 10 < 7 Li/H < 3.5 x 10~ 10 . (5) 

This primordial 4 He abundance constraint includes a statistical uncertainty of ±0.003 
and possible systematic errors as much as +0.01/ — 0.005 with central value of 0.234. A 
recent reinvestigation (with new data) of the linear regression method for estimating the 
primordial 4 He abundance has called into question the systematic uncertainties assigned 
to Y p ([Izatov, Thuan fc Lipovetsky 19961) . Our adopted upper limit to Y p of Eq. (||) is 



essentially equal to the limit derived in their study with la statistical error. 

The upper limit to the lithium abundance adopted here includes the systematic increase 
from the model atmospheres of Thorburn (1994) and the possibility of as much as a factor 
of 2 increase due to stellar destruction. This is consistent with the recent observations 
of 6 Li in halo stars flSmith, Lambert fc Missen 1992| |Hobbs fc Thorburn 1994] ) . We note 



that recent discussion of model atmospheres QKurucz 1995|) suggests that as much as an 



order magnitude upward shift in the primordial lithium abundance could be warranted 
due to the tendency of one- dimensional models to underestimate the ionization of lithium. 
Furthermore, a recent determination of the lithium abundance in the globular cluster M92 
having the metallicity [Fe/H] = -2.25 has indicated that at least one star out of seven shows 
[Li] = 12 + log(Li/H) f=a 2.5 ( |Bocsgaard 19965] ). Since the abundance measurement of the 



globular cluster stars is more reliable than that of field stars, this detection along with 
the possible depletion of lithium in stellar atmospheres suggests that a lower limit to the 
primordial abundance is 3.2 x 10~ 10 < 7 Li/H. 

There also remains the question as to why several stars which are in all respects 
similar to the other stars in the Population II 'lithium plateau', are so lithium rich or 
lithium deficient ( Peliyannis et al. 1995] ; Poesgaard 1996a] , 1996b). Until this is clarified, 
it may be premature to assert that the Population II abundance of lithium reflects the 
primordial value. The primordial abundance may instead correspond to the much higher 
value observed in Population I stars which has been depleted down to the Population II 
lithium plateau. The observational evidence (Deliyannis, Pinsonneault Sz Duncan 1993) for 
a ± 25 % dispersion in the Population II lithium plateau is consistent with this hypothesis 
( peliyannis et al. 1993| ; |Charbonnel 1995| (Steigman 1996b|) . Rotational depletion was 



studied in detail by Pinsonneault et al. (1992) who note that the depletion factor could 
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have been as large as 10. Chaboyer and Demarque (1994) also demonstrated that models 
incorporating both rotation and diffusion provide a good match to the observed 7 Li 
depletion with decreasing temperature in Population II stars and their model indicated that 
the initial lithium abundance could have been as high as 7 Li/H = 1.23 ± 0.28 x 10~ 9 . 

A recent study (|Ryan et al. 1996| ), which includes new data on 7 halo dwarfs, fails 
to find evidence of significant depletion through diffusion, although other mechanisms are 
not excluded. For example, stellar wind-driven mass loss could deplete a high primordial 
lithium abundance of down to the Population II value [Eq. (]5|)] in a manner consistent with 
6 Li observations ( |Vauclair fc Charbonnel 1995| ). Furthermore, it could be possible ( |Yoshii, 
Mathews, fc Kajino 1995] ) that some of the 6 Li is the result of more recent accretion of 
interstellar material that could occur as halo stars episodically plunge through the disk. 
Such a process could mask the earlier destruction of lithium. For comparison, therefore, we 
adopt a conservative upper limit on the primordial lithium abundance of 

7 Li/H < 1.5 x 10~ 9 . (6) 



Finally, the primordial abundance of deuterium is even harder to clarify since it is 
easily destroyed in stars (at temperatures exceeding about 6 x 10 5 K). Previously, limits 
on the deuterium (and also the 3 He) abundances have been inferred from their presence 
in presolar material (e.g., [Walker et al. 199~T|) . It is also inferred from the detection in the 



local interstellar medium (ISM) through its ultraviolet absorption lines in stellar spectra 
( IMcCullough 1992| ; [Linsky et al. 1993| , 1995). The limit from ISM data is consistent with 



that from abundances in presolar material. It has been argued that there are no important 
astrophysical sources of deuterium flKpstein et al. 19"7B ) and ongoing observational attempts 



to detect signs of deuterium synthesis in the Galaxy are so far consistent with this hypothesis 
(see [Pasachoff fc Vidal-Madjar 1989j ). If this is indeed so, then the lowest D abundance 



observed today should provide a lower bound to the primordial abundance. Recent precise 
measurements by Linsky et al. (1995, 1993) using the Hubble Space Telescope implies 

D/H > 1.2 x 10~ 5 . (7) 

We adopt this as a lower limit to the primordial deuterium abundance for the purposes of 
exploring the maximal cosmological impact from IBBN. In addition, we consider the two 
possible detections of the deuterium abundance along the line of sight to high red shifted 
quasars, Eqs. @ and @ as possible limits. 

In order to derive a lower limit to Qb h\o, it is useful to consider the sum of deuterium 
plus 3 He. In the context of a closed-box instantaneous recycling approximation, it is 
straightforward ( |01ive et al. 1990]) to show that the sum of primordial deuterium and 3 He 



- 9- 



can be written 

'X c , 



where A & is the fraction of the initial primordial deuterium still present when the solar 
system formed, is the fraction of 3 He that survives incorporation into a single generation 
of stars, t/230 is the presolar value of [D+ 3 He]/H inferred from the gas rich meteorites, and 
X & /X p is the ratio of the presolar hydrogen mass fraction to the primordial value. These 
factors together imply an upper limit ([Walker et al. 1991 ; Copi, Schramm fc Turner 1995 ) 
of 

V23 P < 1.1 x 10" 4 . (9) 



4. CALCULATIONS 

The calculations described here are based upon the coupled diffusion and nucleosynthesis 
code of Mathews et al. (1990), but with a number of nuclear reaction rates updated and 
the numerical diffusion scheme modified to accommodate cylindrical geometry. We also 
have implemented an improved numerical scheme which gives a more accurate description 
of the effects of proton and ion diffusion, and Compton drag at late times. Although our 
approach is not as sophisticated as that of Jedamzik et al. (1994a), it produces essentially 
the same results for the parameters employed here. We have also included all of the new 
nuclear reaction rates summarized in Smith et al. (1993) as well as those given in Thomas 
et al. (1993). We obtain the same result as Smith et al. (1993) using these rates and 
homogeneous conditions in our IBBN model 

Calculations were performed in a cylindrical geometry both with the high density 
regions in the center (condensed cylinders), and with the high density regions in the outer 
zone of computation (cylindrical shells). Similarly, calculations were made in a spherical 
geometry with the high density regions in the center (condensed spheres) and with the high 
density region in the outer zones of computation (spherical shells). 

In the calculations, the fluctuations are resolved into 16 zones of variable width as 
described by Mathews et al. (1990). We assumed three neutrino flavors and an initially 
homogeneous density within the fluctuations. Such fluctuation shapes are the most likely 
to emerge, for example, after neutrino-induced expansion (IJedamzik fc Fuller 1994Q . We 



use a neutron mean life-time of r„ = 887.0 ( Particle Data Group 1994]) . In addition to the 



cosmological parameter, f2& and fluctuation geometry, there remain three parameters to 
specify the baryon inhomogeneity. They are: R, the density contrast between the high and 
low-density regions; f v , the volume fraction of the high-density region; and r, the average 
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separation distance between fluctuations. 



5. RESULTS 

The parameters R and f v were optimized to allow for the highest values for h\ Q 
while still satisfying the light-element abundance constraints. For fluctuations represented 
by condensed spheres, optimum parameters are R ~ 10 6 and f^ 3 ~ 0.5 ( |Mathews, Kajino| 



fc Orito 1996|) . For other fluctuation geometries, we have found that optimum parameters 



arc: 



R ~ 10 6 ; for all fluctuation geometries 



f^ 3 ~ 0.19; for spherical shells 




0.5; for condensed cylinders 

0.15; for cylindrical shells, 



although there is not much sensitivity to R once R > 10 3 . Regarding f v , we have written 
the appropriate length scale of high density regions, i.e. f^ 3 and f^ 2 for the spherical and 
cylindrical fluctuation geometries, respectively. The variable parameters in the calculation 
are then the fluctuation cell radius r, and the total baryon-to-photon ratio r\ (or h\ Q ). 



5.1. Constraints on f2 6 /i§ 



Figures [T], @, |3(a)| , and [4(a)| show contours of allowed parameters in the r versus rj and 
r versus h\ Q plane for the adopted light-element abundance constraints of Eqs. - (0) 
and for a possible Lyman-a D/H of Eqs. (Q) and ([$]), for the condensed sphere, spherical 
shell, condensed cylinder, and cylindrical shell fluctuation geometries, respectively. The 
fluctuation cell radius r is given in units of meters for a comoving length scale fixed at a 
temperature of kT = 1 MeV. Both of the possible 7 Li limits, Eqs. (gj) and @ which we 
have discussed above, are also drawn as indicated. In order to clearly distinguish the two 
abundance constraints, we use the single and double-cross hatches for the regions allowed 
by the adopted lower (Eq. (|^)) and higher (Eq. @) limits to the 7 Li primordial abundance. 

Even in the IBBN scenario, if the low D/H of Eq. (|3]) ( [Buries fc Tytler 1996| ) is adopted 
as primordial, this range for D/H appears to be compatible with the 7 Li abundance only 
when a higher (Population I) primordial 7 Li abundance limit is adopted, except for a very 
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narrow region of 7] ~ 6 x 10~ 10 and r < 10 2 m. This conclusion remains unchanged for any 
other fluctuation geometries. Therefore, the acceptance of the low (Buries & Tytler 1996) 
value of D/H would strongly suggest that significant depletion of 7 Li has occurred. 



In contrast, adoption of the high D/H of Eq. (0) ( |Rugers fc Hogan 1996a|) as primordial 



allows the concordance of all light-elements. The upper limits to r\ and fib h^ are largely 
determined by D and 7 Li. The concordance range for the baryon density is comparable 
to that for HBBN for small separation distance r. However, there exist other regions 
of the parameter space with optimum separation distance, which roughly corresponds 
to the neutron diffusion length during nucleosynthesis ( |Mathews et al. 1990| ), with an 



increased maximum allowable value of the baryonic contribution to the closure density 
to flbh\ Q < 0.05 for the cylindrical geometry, as displayed in Fig. f4(a)| . This is similar 



to the value for spherical shells as shown in Mathews et al. (1996) and also in Fig. |2| in 
the present work. The condensed sphere limits, however, are essentially unchanged from 
those of the HBBN model. If the primordial 7 Li abundance could be as high as the upper 
limit of Li/H < 1.5 x 10~ 9 , the maximum allowable value of the baryonic content for the 
condensed sphere would increase to fib h\ < 0.08, with similar values for the spherical shell 
( Mathews, Kajino fc Orito 1996|) . For both the condensed cylinders and cylindrical shells, 



the upper limits could be as high as fib h\ < 0.1 as shown in Figs. |3(a)| and |4(a)| . These 
higher upper limits relative to those of the HBBN are of interest since they are consistent 
with the inferred baryonic mass in the form of hot X-ray gas ( |White et al. 1993] ; |White and 



Fabian 1995|) in dense galactic clusters. The acceptance of this consistence, as noted above, 
requires the significant stellar depletion of 7 Li. 



In Figures [3(b)| and |4(b)| , we also show contours for the condensed cylinder and 



cylindrical shell geometries, respectively, but this time with the conventional light-element 
constraints of Eqs. ([|), (5), (7), and (9) as indicated. Since the results for the condensed 
sphere and spherical shell geometries with this set of the conventional abundance constraints 
have already been discussed by Mathews et al. (1996), we do not show those contours here. 
The cylindrical shell geometry of the present work gives the highest allowed value of fib h\ Q . 
Figure |4(b)| shows that the upper limits to rj and fib h\ are largely determined by Y p and 
7 Li. The upper limits for a cylindrical shell geometry could be as high as flbh\ Q < 0.13 
with similar results for the spherical shell geometry ( Mathews, Kajino fc Orito 1996 ). A 



high primordial lithium abundance would increase the allowable baryonic content to as 
high as fib hf < 0.2. The reason that shell geometries allow for higher baryon densities we 
attribute to more efficient neutron diffusion which occurs when the surface area to volume 
area is increased. This allows for more initial diffusion to produce deuterium, and more 
efficient back diffusion to avoid over producing 7 Li. 
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5.2. Observational Signature 

The production of beryllium and boron as well as lithium in IBBN models can be 
sensitive to neutron diffusion. Therefore, their predicted abundances are sensitive to not 
only the fluctuation parameter r, R, and /„ but also the fluctuation geometry ( Boyd fc| 
IKa.jino 19891 ; [Malaney fc Fowler 1989] ; |Ka.jino fc Boyd 1990| ; Terasawa & Sato 1990). 



Figures |5] - [7] show the contours of the calculated abundances for lithium, beryllium and 
boron, respectively in the r versus rj (and r versus fib h\ Q ) plane, the shaded region depict 
is allowed values of r and rj from the light element abundance constraints [cf. Fig. [4(b)|| for 
a cylindrical shell fluctuation geometry. The contour patterns of lithium (Fig. [5]) and boron 
(Fig. 0) abundances are very similar, whereas there is no similarity found between lithium 
(Fig. H) and beryllium (Fig. abundances. 

In order to understand the similarities and differences among these three elemental 
abundances, we show in Figs. || and |9| the decompositions of the A = 7 abundance into 
7 Li and 7 Be and the boron abundance into 10 B and n B. These Figures show also the 
dependence of the predicted LiBeB abundances in IBBN on the scale of fluctuations for a 
cylindrical shell geometry with fixed fib h\ Q = 0.1. This value of fib ^lo corresponds to a 
typical value in the allowable range of rj in Fig. |4(b)| , which optimizes the light element 



abundance constraints, even satisfying the lower 7 Li abundance limit of Eq. ([]). The 
fluctuation parameters f v and R are the same as in Fig. |TF5T. Once the baryonic content 



fib is fixed, the only variable parameter is the separation distance, r. 

As can be seen in Fig. |8|, as the separation r increases, neutron diffusion plays an 
increasingly important role in the production of t and, by the 4 He(t, 7) 7 Li reaction. It 
works maximally around r ~ 10 4 m, which is the typical length scale of neutron diffusion 
at kT = 1 MeV. A similar behavior is observed in the 7 Li(t, n) 9 Be reaction. This reaction 
produces most of the 9 Be in neutron rich environments where t and 7 Li are abundant, as 
was first pointed out by Boyd and Kajino (1989). At other separation distances r in a 
flbh\ = 0.1 model, most of the A = 7 nuclides are created as 7 Be by the 4 He( 3 He, 7) 7 Be 
reaction. In the limit of r = horizon scale, the nucleosynthesis products are approximately 
equal to the sum of those produced in the proton-rich and neutron-rich zones separately 
( |Jcdamzik, Fuller, Mathews fc Kajino 1994|) . The predominant contribution from the 



proton-rich zones makes the 7 Be abundance almost constant at larger r, while both 7 Li and 
9 Be decrease as r increases toward the horizon at any separation distance. 

Figure ^| shows that 11 B is a predominant component of the total boron abundance 
at any separation distance. This is true for almost all values flbh\ Q . It has been pointed 



out flMalaney fc Fowler 1988| ; Applcgatc, Hogan fc Scherrer 1988| ; |Kajino fc Boyd 1990|) 



that most n B is produced by the 7 Li(n, 7) 8 Li(a, n) n B reaction sequence in neutron-rich 
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environments at relatively early times when most of the other heavier nuclides are made. 
Recent measurements of the previously unmeasured 7 Li(a,n) n B reaction cross section ( |Boyd 



et al. 199"2"| ; |Gu et al. 1995| ; |Boyd, Paradellis fc Rolfs 1996] ) at the energies of cosmological 
interest have removed the significant ambiguity in the calculated n B abundance due to 
this reaction. The factor of two discrepancy among several different measurements of the 
reaction cross section for 7 Li(n,7) 8 Li was also resolved by the new measurement ( [Nagai et 



al. 199 1|) . The 7 Li(a,7) n B reaction also makes an appreciable but weaker contribution to 
the production of 11 B in the neutron- rich environment. In the proton-rich environment, 
on the other hand, the 7 Be(a,7) n C reaction contributes largely to the production of n C 
which beta decays to n B in 20.39 min. These facts explain why the contour patterns of the 
lithium and boron abundances in Figs. || and |7| look very similar. 

It is conventional in the literature to quote the beryllium and boron abundance 
relative to H = 10 2 . Hence, one defines the quantity [X] = 12 + log(X/H). In cylindrical 
shell fluctuation geometry the beryllium abundance can take the value of [Be] ~ —3 
while still satisfying all of the light-element abundance constraints and the Population II 
lithium abundance constraint (Figs. |^ and |6]). This abundance is higher by three orders 
magnitude than that produced in the HBBN model with conventional light-element 
abundance constraints. This result is contrary to a recent result with the condensed sphere 
geometry and for a more restricted parameter space ([hhomas et al. 1993 ). Recent beryllium 



observation of Population II stars ( [Rebolo et al. 1988| ; |Ryan et al. 1990| , 1992; |Ryan 199"6] : 



Gilmore et al. 19924 1992b; [Boesgaard fc King 1993| ; Boesgaard 1994| , 1996a,b) have placed 
the upper limit on the primordial 9 Be abundance to [Be] ~ —2, one order magnitude greater 
than the beryllium abundance in the IBBN cylindrical model. 

The calculated boron abundance at the optimum separation distance is essentially 
equal to the value of the HBBN model. However, a high primordial lithium abundance 
would increase the upper limit to h\ Q . In this case, the boron abundance could be one or 
two orders magnitude larger than that of the HBBN model (Fig. [T]). 



CONCLUSIONS 



We have reinvestigated the upper limit to rj and fib h\ Q in inhomogeneous primordial 
nucleosynthesis models. We have considered effects of various geometries. In particular, 
for the first time we consider cylindrical geometry. We have also incorporated recently 
revised light-element abundance constraints including implications of the possible detection 
( [Songaila et al. 1994| ; |Garswell, Science., et, al. 1994 fjHI [fytler fc Fan 1994 ; I'lytler, Fan 
|fe Burlers 1996| ; |Rugers fc Hogan 19964 1996b; |Wampler et al. 1996|) of a high deuterium 
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abundance in Lyman-a absorption systems. We have shown that with low primordial 
deuterium ( [lytler fc Fan 1994| ; |Tytler, Fan, fc Burlers 1996| ), significant depletion of 7 Li is 
required to obtain concordance between predicted light-element abundance of any model of 
BBN and the observationally inferred primordial abundance. If high primordial deuterium 
( IRugers fc Hogan 1996a|) is adopted (Eq. (H)), there is a concordance range which is largely 
determined by D/H, and the upper limit to h\ Q is 0.05. However, with the presently 
adopted (Eqs. (f|), (H), (0), (fD) light-element abundance constraints flSchramm fc Mathews 



|1995| ; |Copi, Schramm fc Turner 1995| ; |01ive fc Scully 19961) , values of Qb h\ Q as large as 0.2 



are possible in IBBN models with cylindrical-shell fluctuation geometry. 

We have also found that significant beryllium and boron production is possible in 
IBBN models without violating the light element abundance constraints. The search for 
the primordial abundance of these elements in low metallicity stars could, therefore, be a 
definitive indicator of the presence or absence of cylindrical baryon inhomogeneities in the 
early universe. 
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Fig. 1. — Contours of allowed values for baryon-to-photon ratio r] (or f^^lo) an d fluctuation 
separation radius r based upon the various light-element abundance constraints as indicated. 
The separation r is given in units of meters comoving at kT = 1 MeV. This calculation is 
based upon baryon density fluctuations represented by condensed spheres. The cross hatched 
region is allowed by the adopted primordial abundance limits with high (Eq. @) and low 
(Eq. (|3])) deuterium abundance in Lyman limit systems and also a higher extreme 7 Li upper 
limit (Eq. @). The single hatched region depicts the allowed parameters for lower 7 Li 
(Eq. (|D) constraint. Note that the 7 Li abundance is the sum of 7 Li and 7 Be. 

Fig. 2. — Same as Fig. 1, but for fluctuations represented by spherical shells. 

Fig. 3(a). — Same as Fig. 1, but for fluctuations represented by condensed cylinders. Adopted 
primordial deuterium abundance constraints are inferred from observations of Lyman limit 
systems (Eqs. (Q) and (3)). 

Fig. 3(b). — Same as Fig. 1, but for fluctuations represented by condensed cylinders. 
Adopted primordial deuterium and 3 He abundance constraints are inferred from observations 
of ISM (Eqs. © and (9)). 

Fig. 4(a). — Same as Fig. 3(a), but for fluctuations represented by cylindrical shells. 

Fig. 4(b). — Same as Fig. 3(b), but for fluctuations represented by cylindrical shells. 

Fig. 5. — Contours of the predicted abundance of lithium (the sum of 7 Li and 7 Be), for 
baryon-to-photon ratio rj (or Qbhl ) and fluctuation separation radius r, in the cylindrical 
shell fluctuation geometry. The shaded region displays the allowed rj—r region from Fig. 4(b). 

Fig. 6. — Same as Fig. 5, but for beryllium, 9 Be. 

Fig. 7. — Same as Fig. 5, but for boron, 10 B + n B. 

Fig. 8. — Lithium and beryllium abundances as function of proper separation distance in 
units of meters comoving at kT = 1 MeV for fixed Qbh\ Q = 0.1. Refer to the abundance 
scales in l.h.s. for lithium and r.h.s. for beryllium. 

Fig. 9. — Same as Fig. 8, but for boron. 



